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AN ATLAS FOR INTERPRETING 7-RAY PULSAR LIGHT CURVES 
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ABSTRACT 

We have simulated a population of young spin-powered pulsars and computed the beaming pattern and light 
curves for the three main geometrical models: polar cap emission, two-pole caustic ("slot gap") emission and 
outer magnetosphere emission. The light curve shapes depend sensitively on the magnetic inclination a and 
viewing angle £. We present the results as maps of observables such as peak multiplicity and 7-ray peak 
separation in the (a, £) plane. These diagrams can be used to locate allowed regions for radio-loud and radio- 
quiet pulsars and to convert observed fluxes to true all-sky emission. 
Subject headings: gamma-rays - pulsars - stars: neutron 



1. INTRODUCTION 

With the successful launch of the Fermi Gamma-ray Space 
Telescope, formerly GLAST, many pulsars are expected to be 
detected in the 7-ray band. This is an important opportunity 
since ~GeV 7-ray emission dominates the observed electro- 
magnetic output of young pulsars. Also, some pulsars are 
expected to be detected in the 7-ray band but undetectably 
faint in the radio. The natural interpretation is that the radio 
and 7-ray antenna patterns differ. Indeed, the observed pul- 
sars have very different radio and 7-ray pulse profiles. Since 
the observed pulse profile is simply a cut through this antenna 
pattern for the Earth line-of-sight angle £g, interpretation of 
the radio and 7-ray profile requires a beaming model. When 
the spin geometry is known, observed pulses can be used to 
select the best model. Conversely, for a given model, the pulse 
pattern can restrict the spin orientation. Finally, beaming cor- 
rections can be used to relate the observed line-of-sight flux 
to the full-sky emission. Thus, aided by beaming models, 7- 
ray pulsar data can give a new window into the energetics and 
evolution of neutron stars in the Galaxy. 

There are three main geometrical models for the 7-ray pulse 
beaming, each tied to locations in the magnetosphere where 
force-free perfect MHD (E ■ B = 0) conditions break down 
and the uncanceled rotation-induced EMF E w ril x B/c 
causes particle acceleration and radiation. These are the so- 
called magnetospheric gaps. The first 7-ray (and radio) pul- 
sar models traced emission to acceleration zones on the field 
lines a t very low altitudes above each magnetic pole (al- 
though Daughe rtv & Hard ing 1996 argued that such "polar 
cap" model emission could extend to ~ 2i?»). Another major 
class of models posits "Holloway" gaps above the null charge 
surface, ext ending toward the light cy linder (the "outer gap " 
model; see ICheng. Ho. & Rudermanlll986t lRomanilll996l) . 
More recently, it has been argued that t he true null ch arge 
location depends on currents in the gap (lHirotanill2006l) . al- 
lowing the start of the outer gap to migrate inward towards 
the stellar surface. Conversely, it has been argued that accel- 
eration at the rims of the polar caps may extend to very high 
altitude (Muslimov & Harding 2004), effectively pushing the 
polar cap activity outward. Together these effects suggest that 
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emission can occur over a large fraction of the boundary of the 
open zone. One model that has been suggested for this geom- 
etry makes low-altitude emission visible from one hemisphere 
and higher-altitude emission (i.e., above the null charge sur- 
face) visible from the other. This is the "two-pole caustic" 
(TPC) model dDvks & Rudakll2003l) . which is thus intermedi- 
ate between the polar cap (PC) and outer gap (OG) pictures. 

2. SIMULATION METHOD 

We start by building a library of pulsar antenna patterns. 
We follow stand ard practice in basing our beaming on the re- 
tarded potential (Deutsch 1955) magnetic dipole fields. This 
has the appreciable advantage of analytic simplicity, allowing 
rapid computation s. Numerical mo dels of force-free magne- 
tospheres do exist (ISpitkovsky 2006), but such models do not 
have the resolution required for the sort of beaming compu- 
tations pursued here and in any case do not, by definition, 
include acceleration gaps. Still, it should be remembered that 
magnetospheric charges and currents may well have an ob- 
servable effect on the field structure and pulse shapes com- 
puted here. 

We model the last closed field line surface (those field lines 
tangential to the speed of light cylinder at r± — Tlc = 
cP/2ir), tracing these field lines to the polar cap surface at 
the stellar radius r. t . The relevant magnetospheric features in 
the OG and TPC models are tied to the light cylinder and lo- 
cations identified as fractions thereof; thus, the period of the 
pulsar is an unnecessary parameter in mapping the structure 
of the magnetosphere. Models are computed across all possi- 
ble magnetic inclinations a. We follow the field line structure 
to <10~ 4 rLcs allowing precision mapping of the PC emis- 
sion region even for pulsar periods >1 s. A single magneto- 
sphere model can serve a range of pulsar periods by truncating 
the grid at larger fractions of r^c for shorter period pulsars, 
thereby maintaining the same physical distances. General rel- 
ativistic effects make small changes for r < 3r* but are not 
included in the present sums. 

In all models, the acceleration gap is inferred to arise near 
the open zone boundary tied to the last closed field lines. For 
the PC picture, we simply follow radiation emitted tangen- 
tially to these field lines at al titudes <IR*. In th e original 
definition of the TPC model dDvks & Rudakll2003l) the emis- 
sion was similarly placed on the last closed field lines but ex- 
tended to an altitude r < r^c and a perpendicular distance 
r± < 0.75rj,c> where r^c is the perpendicular distance to 
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FIG. 1. — Sample antenna intensity patterns for a = 65° with the horizontal axis of pulsar phase (<j>) and vertical axis of viewing angle (Q. The TPC model 
for w = 0.05 is shown in green, the OG model for w = 0.1 is shown in red and the PC emission site is shown in blue. The cyan lines show the locus of the 
possible high altitude (r = 500 km, here for P = 0.2 s) radio emission, with the radiating front half shown solid and the back half dashed (see §3.3). 



the light cylinder from the rotation axis. In the outer gap (OG) 
model the emission is from a full set of open field lines asso- 
ciated with the closed zone surface, but pair production and 
radiation are expected to start above the null-charge surface 

[fi ■ B{R NC ) = 0]. 

In any physical model, however, there should be an open 
zone vacuum region separating the surface of last closed field 
lines from the pair formation front, whic h defines the rad iat- 
ing surface of the gap. The OG model of Roman| d 19961) de- 
fined gap thickness w as the fraction of the angle from the last 
closed field line to the magnetic axis that remains a charge- 
starved vacuum with a large acceleration field. This char- 
acteristic gap width w is proportional to the 7-ray efficiency 
(see below). A physical TPC model sho uld also have a finite 
gap th ickness; the "slot gap" model of Mus limov & H arding 
(2003) provides a possible physical foundation for the TPC 
picture. These authors give an estimate for the dependence 
of gap thickness on spin period P and surface magnetic field 
B; again the thickness increases for pulsars putting a larger 
fraction of the spindown power into pulsed 7-rays. Thus, for 
both models, we assume that the gap thickness can be related 
to a heuristic 7-ray luminosity 

L 7 « V E SD « C x ( 10 33^ s _ 1 ) 1/2 x 10 33 erg (1) 

with C, a slowly varying function of order unity which must 
come from a detailed physical model. This law, which as- 

* — 1/2 

sumes a 7-ray efficiency [w =)rj oc E SD , is natural in mod- 
els that maintain a fixed voltage drop across the acceleration 
gap. Clearly, as w — ► 1 this simple geometrical approxima- 
tion should break down. The original TPC formulation has 
w = 0. In a full physical model, the radiating pair formation 
front would have a finite thickness and the radiation pattern 
of the particles would have a finite width; these effects are not 
treated here but would serve to smooth out the beams com- 
puted in the basic geometrical model. Sample antenna pat- 
terns (sky maps of the pulsar beam) are shown for the three 
models in Figure Q] 

In each case, the antenna pattern is built up by Monte Carlo 
simulation. Radiation is emitted uniformly along the gap sur- 
face. The photons are aberrated, the time of flight to a distant 
observer is computed and the emission is assigned to a bin 
of p ulsar phase and co-latitude (th e sky map antenna pattern; 
see iRomani & Y adigarogiul 119951: iDvks. Harding. & Rudakl 
l200lT We normalize this antenna pattern by summing emis- 



sion at all pulsar phases <p and along all lines-of-sight £ and 
scaling by the heuristic 7-ray luminosity (Equation 1). 

Although we focus here on the beaming characteristics 
rather than the population, to obtain a realistic sample of pul- 
sar spin properties we simulate a Gala ctic population of young 
pulsars following the method of IWeltevrede & Johnston 
(2008a]). P ulsars are assigned a P and P fr om the parent dis- 
tribution of Weltevrede & Johnston] d2008al) . Since we are in- 
terested in energetic pulsars, we retain objects with E > 10 34 
(i.e., a max characteristic age t = P/2P w 10 6 y), drawing 
until we have a sample of ^175,000 pulsars. This is ~12x 
the true Galactic density (given a birthrate of ~1.5 per cen- 
tury); the over-density serves to decrease statistical fluctua- 
tions. Magnetic inclinations a and Earth viewing angles £ 
were drawn isotropically. Given the low upper limit on age, 
we ignore any possible magnetic alignment, which, if present, 
appears to occur on longer time scales. For each drawn pulsar 
the appropriate 7-ray light curve is computed for the drawn 
line-of-sight ( using the appropriate antenna pattern from the 
library of sky maps. 

While we do use a realistic pulsar population, we do not 
report here on predictions for the statistical distributions of 7- 
ray pulse shape, luminosity and Galactic location. These, and 
predictions for the relative numbers of radio-loud and radio- 
quiet pulsars expected in the LAT sky survey, are deferred to 
a subsequent paper. We also do not treat here the millisecond 
pulsar population, as with small light cylinder radii the differ- 
ence between inner magnetosphere and outer magnetosphere 
emission zones is less stark. 

2.1. Radio Emission 

Radio beam models have a long history but are, if anything, 
even more heuristic than those for the GeV 7-rays. Here we 
consider two scenarios. One is simple low-altitude emission, 
which assumes a circular cone of radio emission along the last 
closed field lines with an angular width 

2p = 10.8P~ 1/2 degrees (2) 

centered on the magnetic axis dRankinl [19931 ; [Gouldl 1 1 9941) . 
Here, p is the half-opening angle of the radio cone; the scaling 
is well supported by observations, which also agree with the 
beam being roughly circular, as expected from dipole field 
line structure at moderate altitudes. To model this scenario 
in our baseline computations we assume radio observability 
whenever the line of sight passes within p of the magnetic 
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FIG. 2. — An atlas of pulse profile properties for three 7-ray pulsar models. Top row: pulse complexity (peak number); broad peaks lie in zone enclosed by the 
cyan lines. Middle row: pulse width A, the maximum separation of the major peaks. Bottom row: the "flux conversion factor" /o for correction to all-sky flux. 
For the PC model, the small ranges of A and /n require re-normalization of the color scale for those panels. 
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axis. The radio pulse peak is assumed to be at the same phase 
as this axis. 

Recent evidence, however, suggests that young, Vela-like 
pulsars have a fixed emission altitude of ~ 5 00 — 1000 km or 
~ 30— 70-fL dKarastergiou & Johns ton 2007). Since these are 
the dominant denizens of the 7-ray pulsar population, we also 
(see §3.3) consider the effect of such wide beams at higher 
altitude. This wide beam emission is assumed to be patchy, 
with some lines of sight crossing the radi o zone but receiving 
negligible flux dLvne & Manch ester 1988). This is difficult to 
treat, except statistically, and is only discussed qualitatively in 
this work. As we shall see in §3.3, it is useful to consider the 
possibility that the leading portion of the radio zone has higher 
occupancy than the trailing half. We thus also compute radio 
pulse observability for emission only along the leading edge 
of a 500 km altitude radio emission zone (Figure [TJ. Aberra- 
tion and light travel effects are computed for all radio pulses 
as well. 

2.2. Pulse Profile Characterization 

We wish to summarize the properties of the 7-ray pulses 
of the simulated population. We start by flagging the major 
peaks for each light curve, generated from cuts at a given £ 
across antenna patterns such as those in Figure 1 . Examples 
of light curves for these model s can be seen in Figure 5 of 
iDvks. Harding. & Rudak ( 2004) for the particular case a = 
60°. In the Appendix we show two summary figures of the 
light curves for a range of a and £ (and w) for the OG and 
TPC models. 

We have created a peak-finding routine that runs on the 7- 
ray light curves produced from the models. The light curves 
are searched for maxima that rise above the local minimum 
by least 20% of the global light curve maximum. Most non- 
thermal pulsar peaks are relatively narrow, so we flag "sharp" 
peaks by requiring that the the full width at half maximum 
(FWHM) of the qualifying peak (i.e. measured at 10% of the 
full pulse amplitude below the peak maximum) have A(f> < 
0.1. The bright EGRET pulsars with high S/N light curves 
(Crab, Vela and Geminga) all show two peaks using this cut. 
However, wider maxima may also be of interest, so we addi- 
tionally flag "broad peaks" with FWHM 0.1 < A<p < 0.3. 
These tend to appear when either ( or a is small. Note that 
some pulsars with both small £ and a have appreciable 7- 
ray flux but show only shallow sinusoidal variations with no 
strong "peaks" as defined by the above algorithm. 

Two additional key values regarding peak locations are also 
recorded. First, we record the separation from the first peak 
to the last peak. Additionally, we measure the lag of the first 
peak from the radio pulse. To help the reader understand 
the results of this peak flagging, in the Appendix figures we 
record the number of peaks, the broad peaks and the pulse 
width (in %) as determined by this routine for each of the 
plotted light curves. 

Finally, we compute fn, the ratio between the simulated to- 
tal emission of the pulsar over all observer angles ( and the 
total emission if the observed flux at were emitted uni- 
formly across the sky. 

3. RESULTS 

We have collected the geometrical pulse properties of the 
three models into an "atlas" (Figure |2} that compactly sum- 
marizes many complex pulse profiles. We plot simulated pul- 
sars in the (a, Q plane, color coding the points to show the 
model properties. Small dots show objects detected in 7-rays 



but not in the radio. Radio detectable objects tend to lie at 
small (3 — C, — a, the magnetic pole impact angle. The dashed 
lines show the low-altitude standard radio detection window 
for P — 0.2 s, the median for the observable 7-ray pulsar 
population. Of course, these are not strict boundaries; pulsars 
with large P have smaller radio pulse widths and may remain 
radio invisible at smaller j3 while especially fast pulsars may 
be radio detected well beyond the plotted lines. There are 
radio-detectable, 7-invisible objects as well. These are gener- 
ally found at small a, £ and are not plotted here. 

The top row shows the pulse complexity, the number of ma- 
jor pulse peaks, with broad peaks appearing in the cyan curve. 
The second row shows, for each model, the phase interval A 
between the most widely separated major peaks. These plots 
may be used to predict the 7-ray properties for a given model 
and known orientation angles. Conversely, for an observed 
7-ray pulse profile one may restrict the allowed viewing an- 
gles for a given model or even rule a model out completely. 
Of course, objects with well-constrained viewing geometry 
and observed profiles are particularly useful as they allow the 
most strict discrimination between pulse models. More gen- 
erally, the distribution of observed pulse shapes is different 
between the models. Once one assumes an underlying pulsar 
distribution, such population analysis is also a very powerful 
discriminant (Watters et al. 2009). 

The final row of the atlas depicts a flux correction factor fn, 
useful in calculation of 7-ray efficiencies. The pulse profile 
observed along the line of sight at £ for a pulsar with mag- 
netic inclination a is F 1 (a 1 £, 0), which we can compute for 
a given pulsar model. This means that the observed phase- 
averaged flux F bs for the Earth line-of-sight C,e is not neces- 
sarily representative of the flux averaged over the sky. This 
has al so been emphasized by lHarding. Grenier & Gonthien 
(2007). Here, we can use our model to compute the required 
correction fn, so that the true luminosity is 



L 7 = 4irf n F ohs D 2 . 
where D is the distance to the pulsar and 

, f ( , v //F 7 (a, C,0)sin(OdC<ty 



(3) 



(4) 

2jF 1 (a,C E ,W ' W 

It has become traditional to assume that the 7-ray beam cov- 
ers an area of 1 sr uniformly, so that fn — l/(47r) « 0.08; 
however, more modern pulsar beaming models have fn ~ 1 
for many viewing angles. 

The flux correction factor is crucial for estimating the 7-ray 
efficiency 

r\ = L-f/Esn oc f n . (5) 

Thus, it appears that many recent papers have underestimated 
the true 7-ray efficiency (for modern pulsar models) by an 
order of magnitude. 
From our heuristic luminosity law, we expect that the effi- 

• —1/2 

ciency evolves as 77 oc E SB . In general, one would expect a 
larger fraction of the gap to remain charged-starved for more 
efficient pulsars, with w oc r\. In figures [3] and [4] we see how 
the pulse width and sky covering factor evolve with w = ?/. In 
particular, for low-Fgo, large- w, high-efficiency pulsars, we 
see that fn decreases. Thus, the correction of observed phase 
average flux to true sky flux becomes a function of spindown 
power. As w saturates at a large value of order unity, which 
for this calculation is near Fgo = 10 33 erg s _1 , no radiating 
surfaces bound the gaps and 7-ray production ceases. The 
exact value depends on magnetic inclination a. 
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Main peak separation 

•0 - 0.1 
•0.1 - 0.2 
• 0.2 - 0.3 




Magnetic Inclination 



FIG. 3. — Beam shape evolution with efficiency for the TPC model. The three rows show peak complexity, peak separation and flux conversion factor, as in 
Figure[2] The gap width (and hence 7-ray efficiency) increases left to right. 
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For the convenience of other researchers, we have collected 
approximate fitting formulae for fn, the flux correction factor. 
Because of the complexity of the beam patterns, it is always 
best to use the atlas figures to determine ffi for a pulsar with 
known parameters. However, statistical statements about lu- 
minosity and rough population sums can benefit from analyt- 
ical approximations. For both the TPC and OG models there 
is a change in the beaming behavior at large a and (. In this 
region one sees emission from the hollow cone above the null 
charge surface that dominates the pulsation for the classical 
OG and often provides the late pulses in the TPC picture. We 
call this "case I" and thus define an approximate boundary for 
such objects as 

C > Ci = (75 + 100w) - (60 + l/w)(a/90) 2(1 ~ tu) , (6) 

with ( in degrees. Any pulsars with £ < (j compose case 
II. For the OG scenario, < 10% of the modeled pulsars have 
any emission at smaller £, so case II is not a major contributor 
to the population. Such objects do, however, exist: princi- 
pally pulsars with small w and large inclination a that are 
observed at small (. The emission comes from diverging 
high-altitude field lines producing faint 7-ray emission (of- 
ten with no strong pulse) over much of the sky. In the TPC 
model, however, there is substantially more emission in case 
II, since the flux at these angles arises from low-altitude field 
lines viewed relatively near the polar cap. These light curves 
tend to show shallow sinusoidal pulses. They are also missing 
the caustic peaks (sharp pulses) dominating the case I objects. 
Thus, in both the TPC and OG models, Equation 6 defines 
a boundary past which the pulse emission is fainter and /□ 
is larger. Figures |2][5] and the light curves in the Appendix 
(Figures l9land[T0b make these patterns clear. 
For the TPC model, we find in the case I zone 

/ Q » 0.8 + 1.2(0.3 -«j)cos(2/3), (7) 

where (3 = Q — a. This applies for pulsars with w > 0.05 
and gives estimates accurate to 10% for 76% of the pulsars 
and to 20% for 90% of the pulsars. In the case II zone the flux 
correction factors are generally larger, with 

fa « 0.3 + 1.5(1 - «>)[1 + (C - C/)/90]. (8) 

This region also has more complex behavior, so it is better 
to use the atlas figures if possible. The overall accuracy for 
the TPC model estimates (case I and case II zones together) 
is better than 10% for 62% of the modeled objects and within 
20% for 82% of all modeled pulsars. 
For the OG we consider the case I zone only, finding 

fn ~ 0.17 - 0.69u> + (1.15 - 1.05w)(a/90) 19 (9) 

for £ > 60°. A ^-dependent term appears for pulsars viewed 
closer to the spin axis, so that 

fn w 0.17 - 0.69u; + (1.15 - 1.05w)(a/90) 1 ' 9 

-1.35(2/3 -C/90) (10) 

when ( < 60°. These give approximations accurate to better 
than 10% for 70% of the simulated case I pulsars with w > 
0.05 and are better than 20% for 92% of such pulsars. The 
atlas figures give a guide to regions where exceptions occur. 

3.1. Determining Individual Pulsar Geometries 

Many young energetic pulsars produce bright, resolvable 
X-ray pulsar wind nebulae (PWNe), the beam dump of rela- 
tivistic particles and fields flowing out of the magnetosphere. 



Often, high-resolution X-ray images of these PWNe can be in- 
terpreted in terms of cylindrically symmetric structures, with 
a mildly relativistic (v « 0.3 — 0.7c) bulk outflow Doppler- 
boosting the intensity along the radial flow. The inference is 
that the symmetry axis traces the pulsar spin, and fits to such 
models dNg & Romanil l2008h can provide relatively model- 
free estimates of Thus, for such pulsars, the family of 
allowed models is given by horizontal slices across the atlas 
panels. 

Since we are concerned with detected 7-ray pulsars, one 
will inevitably have an estimate of the pulse width A. This is, 
then, the most robust observable, and provides a key to deter- 
mining the rest of the model parameters. Thus, the prescrip- 
tion for solving a pulsar for a given model is to i) use Esb 
(Equation 1) to estimate w; ii) go to the atlas pulse width plot 
for the closest w and locate regions showing a suitable A; Hi) 
use an external value for £e, when available, to locate the rel- 
evant point in the A band and confirm that that the pulsar is 
radio loud or quiet, as appropriate; iv) use the £ and a values 
for this point to read off other pulsar properties (Figures 2-4). 
Since the A plot is the "key" to the atlas, and since there are 
large regions at high a and £ with A ~ 0.4 — 0.5, especially 
for the TPC model, we provide in Figure 5 an expansion of 
the lower right quadrant of the evolving A panels with finer 
color bands. This should be used to determine £ and a when 
possible, which can be transferred to Figures 3 and 4. 

There is also a well-developed phenomenology connecting 
the sweep of linear polarization within the radio pulses to cuts 
across a simple magnetic pole. In the classic rotating-vector 
model (RVM, Radhakris hnan & Cookdfl969l) . the emission 
is assumed to have a polarization vector locked to the dipole 
magnetic field (projected on the plane of the sky) in the emit- 
ting region. The shape of the polarization sweep depends 
most sensitively on the value sin(ce)/sin( / 9). As (3 generally 
takes on small values (or else the object would not be detected 
in the radio in the first place), this method tends to give poor 
constraints on a and tighter constraints on j3. Thus, in prin- 
ciple, X-ray and radio data together can provide both £ and a 
and therefore a prediction of the 7-ray pulse shape. 

In practice, the picture is often not this simple. One pos- 
sible complication arises if the altitude of the radio emission 
is large, where sweep-back and tim e-delay can alter the naiv e 
rotating-vector model predictions (Everett & Weisberg 2001). 
Finally, we should reiterate that our models are purely geo- 
metrical (albeit motivated by physics) and locked to the vac- 
uum dipole field. Field-line perturbations, currents and vari- 
ation in emissivity along the gap can all perturb the profiles. 
Nevertheless, this atlas does capture the key differences be- 
tween the models and provides a basis for broad-brush inter- 
pretation of newly discovered pulsars. 

3.2. The Polar Cap Model 

The polar cap picture naturally predicts that only pulsars 
seen at small (3 will be seen in the 7-rays. With the radio 
emission originating above the surface 7-rays (which are at 
r < 2r*), one would always expect radio emission bracket- 
ing the GeV pulse. "Radio-quiet" objects, such as Geminga, 
are not possible in this picture unless the radio cone is very 
patchy and emission is completely absent along Qe- "Inter- 
pulsars," with emission visible from both radio poles, are ex- 
pected at the lower right (a ~ C ~ 90°), where the left panel 
of Figure 2 shows triple- and quadruple-peaked light curves 
for this model. Note that there is also a small collection of 
triply-peaked light curves at a ~ ( ~ 45° stemming from 
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4. — Outer Gap model beam shape evolution as a function of gap width, as in Figure [5] 
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FIG. 5. — The pulse width for TPC and OG models at several w. This figure should be used when possible to locate accurately the a and f for a given pulsar, 
for reference back to Figures[5]and|4] 

the notch structure observe d in the outline of the polar cap f or 
inclination angles a — 45° dDvks. Harding. & Ru dak 2004). 
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FIG. 6. — Pulsar 7-ray peak separation vs. radio lag. Green points are 
radio-quiet 7-ray pulsars, blue points are radio-loud 7-ray pulsars. 



The wide double pulses, with A = 0.35 — 0.45, seen 
for the bright EGRET pulsars pose another challenge to the 
PC model. While a double-pulse light curve is the domi- 
nant mode in the PC zone, a wide separation can only be 
obtained for very small a and a small range of Q. Unless 
pulsars are born highly aligned, this is statistically very im- 



0.2 0.3 0.1 

Radio Lag (6) 

FIG. 7. — As Figure[6] but assuming high-altitude radio emission from the 
leading half of the radio zone. 

probable. Additionally, model-independent measurements of 
£ made through the analysis of PWN tori mentioned above 
have consistently found values £ > 40°, well beyond the 
range for which the polar cap model can produce large peak 
separations. 

These geometrical arguments, as well as considerations of 
flux and spectrum, imply that the bulk of the observed 7-ray 
pulsar emission seen by EGRET could not have come from a 
low-altitude cap. However, if a surface gap does have active 
pair production, some 7-rays are nevertheless expected from 
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this region; very sensitive Fermi LAT observations may be 
able to uncover such a component. 

3.3. Observed Pulsars 

The atlas lets one read off the expected 7-ray pulse prop- 
erties for known a and £ or estimate from A values for un- 
determined a and £. It is also interesting to compare the 
observed properties of specific pulsars with the models' pre- 
dictions. Table Q] lists the observed properties of the well 
known EGRET pulsars and two ne w detections, J202 1+3651 
from AGIL E (iHalpem et alJ 120081) and CTA 1 from Fermi 
(Abd oetal.1 12008)" Column 4 lists the externally derived val- 
ues for the line-of-sight angle £, when available. 

Note that for PSR B 1055 -52 the radio-derived C accords 
well with the value predicted by the 7-ray pulse width A in 
the outer gap picture. The TPC picture does have a solution 
matching the observed A (albeit only for w much lower than 
inferred from Esb), but £ does not agree with the radio value, 
and additionally a is too small to produce the observed radio 
interpulse. For the two radio-quiet pulsars Geminga and PSR 
J0007+7303 in CTA 1 we do not have C estimates. For CTA 
1, both OG and TPC models can match the observed A. For 
Geminga, with A = 0.5, the TPC model has solutions along 
both the £ and a axes. For the OG model, Geminga-type A 
are only natural at small a and large £, where one will almost 
certainly miss the radio beam. 

In addition to the 7-ray pulse width A, for radio-detected 
objects, the lag S of the first 7-ray pulse from the radio peak 
is a second convenient observable. In Figure [6] we have plot- 
ted in blue a sample of simulated 7-ray pulsars, in the A — S 
plane, assuming that the radio peak is centered along the mag- 
netic axis. If the viewing angle £ is too far from the pole (i.e., 
j3 large) then no radio pulse will be detected. In these cases we 
have plotted the pulsars in green, with 5 the phase lag from the 
closest approach to the radio pole. For the TPC model, many 
of these pulsars have a first 7-ray peak leading the radio pole. 
Some known 7-ray pulsars with well defined pulse profiles 
have been placed on the figure. For the two "radio-quiet" pul- 
sars Geminga and PSR J0007+7302, no 5 is available, so the 
allowed regions are shown as horizontal bands. PSR B1509- 
58 was seen at MeV energies by CGRO, but was not convinc- 
ingly seen by EGRET; its broad peak has a small but poorly 
defined A, so it appears as a diagonal band. 

It is clear from Figure[6]that many of the observed radios- 
ray pulsars have a larger radio lag than expected in either the 
TPC or OG model if the radio pulse is centered on the mag- 
netic axis. This suggests that some ingredient is missing from 
our basic geometrical model. One possibility is that magne- 
tospheric currents sweep the 7-ray emitting field lines back 
to later phase. Alternatively, an extra lag may be connected 
with th e posited high altitude radio em ission of the young 
pulsars dWeltevrede & Johnsto n 2008b). If we invoke patchy 
emission and make the additional (new) assumption that the 
patches are most active on the leading half of the radio zone, 
then aberration and field-line spreading on this leading half of 
the cone can drive the high-altitude radio peak to much earlier 
phases. The range of allowed lags is shown in Figure [7]for a 
fixed altitude of 500 km. Note that the phase lead i5 then de- 
pends on the pulsar spin period, since this fixed altitude is a 
variable fraction of tlc- These are the model tracks in Figure 
|7] Large (3 will result in a smaller increase in S as the ra- 
dio cone is cut with a small chord, further from the magnetic 
axis. Again for the green (radio-quiet pulsar) dots, we plot the 
S from the closest approach to the radio pole. As expected, the 



wider radio pulse allows fewer green (radio non-detected) ob- 
jects. 

Finally, we should compare the model pulsar efficiencies 77 
with the observed 7-ray fluxes, using our computed /q. Of 
course, this is a purely geometrical comparison, ignoring real 
differences in 7-ray spectrum and the radiation processes tap- 
ping the primary electron energy. Nevertheless, the differ- 
ences in ffi are substantial (and often completely ignored in 
discussions of 77). For gaps with near-constant brightness and 
pulsars with well determined distances, the fn variations can 
be the biggest unknown, making our treatment useful. 

Table [TJ contains the estimated /q ranges inferred for the 
TPC and OG models. In Figure [S] we plot the inferred ef- 
ficiencies of the detected 7-ray pulsars, after correction with 
/n to simulated all-sky emission for the two models. The 
distance uncertainties generally dominate the range in the de- 
rived i 7 . The solid line shows complete conversion of spin- 
down power to 7-rays and the dashed line shows the heuristic 
efficiency law assumed above. 

In Figure [8] two objects near log(E) = 36.5 deserve com- 
ment. The first is PSR J202 1+3651 in the "Dragonfly" neb- 
ula. This has a very large dispersion measure, and the cor- 
responding 12.4kpc distance would imply > 100% efficient 
7-ray production for both TPC and OG pictures. However, 
IVan Etten. RomanjT& Ng (2008) argue that a variety of pulsar 
and PWN measurements are more consistent with a distance 
of < 3 kpc. In fact, the thermal surface emission has a best fit 
distance of d = 2.1 kpc; d ~ 1.6 kpc would bring this pulsar 
into agreement with the heuristic efficiency law. The origin 
of the large DM for this pulsar is thus very puzzling. The 
next outlier is PSR B1706— 44, at si milar E, where we hav e 
adopted the DM distance of 2.3 kpc (ICordes & Lazioll2002T) . 
Again, smaller distances are implied by fits t o the neutron 
star th ermal emission and the PWN energetics (Roma nTet al.1 
2005), with preferred values at d < 1.8 kpc. However, to 
match the efficiency law would require distances close to 
1 kpc; this is substantially below the statistica l lower bound on 
the di stance estimated from H I absorption dKoribalski et al.1 
1 1995b . Thus, the Earth line-of-sight flux of these two Vela- 
type pulsars seems larger than predicted from their spindown 
luminosity and the beaming model. Improved flux measure- 
ments from the LAT and improved distance estimates (e.g., a 
radio interferometric parallax for PSR B 1706 —44) would be 
particularly valuable in tightening up this argument. 

The most interesting region of Figure |8]is the low-E range, 
where the modeled efficiencies approach unity. For these 
large w pulsars, fa begins to drop significantly below unity, 
especially in the OG model, decreasing the effective L 1 . This 
effect dominates as pulsars approach the 7-ray death zone, 
near the end of their lives as high-energy emitters. The small 
fa values prevent the apparent catastrophe of L 1 > E for old 
pulsars. Conversely, while a small fraction of such pulsars, 
with their 7-beams closely confined to the spin equator, are 
visible along the Earth line of sight, those that are seen have 
relatively high fluxes and are detectable to relatively large dis- 
tances. Geminga is an excellent example. 

4. DISCUSSION & CONCLUSIONS 

Our study allows a quick visual summary of model predic- 
tions, which may be compared with pulsar discoveries antici- 
pated from the LAT. We have already seen that the classic po- 
lar cap (PC) picture is difficult to reconcile with the data. Both 
the TPC and OG pictures, however, seem viable. Certain ob- 
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TABLE 1 

Geometry and Beaming of Bright EGRET Pulsars 



Name 


log(_E S D) 


w 


d (kpc) a 


c b 


A 1 ' 


5 C 


C«TPC d 


(tpc a 


/o(TPC) d 


"OG d 


CoG d 


/ S2 (OG) d 


Crab 


38.7 


0.001 


2.0 ±0.2 


63 


0.40 


0.111 


55-60 


63 


1.1 


70 


63 


1.0 


Vela 


36.8 


0.01 


n 9S7 -0.017 
u ' z °'+0.019 


64 


0.426 


0.117 


62-68 


64 


1.1-1.15 


75 


64 


1.0 


B1951+32 


36.6 


0.02 


2.0 ± 0.6 




0.48 


0.135 


55-90 


55-90 


1.0-1.25 


60-90 


70-85 


0.75-1.1 


B 1706-44 


36.5 


0.02 


2.3 ±0.7 


53 


0.2 


0.25 


45-50 


53 


1.05 


47-49 


53 


0.7-1.0 


J202 1+3651 


36.5 


0.02 


3.0 ± 1.0 


86 


0.465 


0.165 


70 


86 


1.1 


70 


86 


1.05 


CTA 1 


35.7 


0.05 


1.4 ±0.3 




0.20 




55-70 


23-45 


0.8-1.3 


20-50 


60-75 


0.3-0.9 


Geminga 


34.5 


0.18 


n 9 c-0.062 
U - ZJ +0.120 




0.50 




30-80,90 


90,55-80 


0.7-0.9,0.6-0.8 


10-25 


85 


0.1-0.15 


B1055-52 


34.5 


0.18 


0.72 ±0.15 


67 


0.22 


0.32 


60 c 


50 


0.85 


78 


67 


0.55 



Distance references (in order): traditional for Crab; Dodson et al. (2003); Strom & Stappers (2000); Cordes & Lazio (2002); Van Etten. Romani. & Ng 



12008); Abdo et al. (2008); Fahertv, Walter, & Anderson (2007); Cordes & Lazio (2002). b From torus fits to X-ray images (Ng & Romani 2008) when avail- 
ab le; B1055-52fro m radio RVM fit <Lvne & Manchester 1988|). c Derived from light curves in lFierrol fl99S) except J2021+365f (H alpern et al.l2008h and CTA 
1 lAbdo et all200SD . d Estimated from Figures |5]and|4] e Requires w < 0.1 




34 36 38 

log spindown energy (ergs/s) 

FIG. 8. — Derived 7-ray luminosity, corrected (TPC - crosses, OG - cir- 
cles) from the observed phase-averaged flux to all-sky, vs. the spindown 
luminosity. Error bars include both the d and /n uncertainties (Table 1). The 
dotted error bar at large L 7 represents the 7-ray luminosity of J202 1+3651 
assuming the DM distance of 12.4 kpc. 

servations would clearly discriminate. For example, detection 
of 7-ray pulsars with a + £ < 90°, but with nearly constant 
intensity through the period, would give a clear preference for 
the TPC model. In contrast, the prevalence of light curves 
with simple double pulses separated by 0.05 < A < 0.3 is 
a hallmark of high-altitude OG emission. The distributions 
in peak widths and in the relative number of radio-loud and 
radio-quiet pulsars are also a strong discriminant; we will ad- 
dress these in a following paper concentrating on the popula- 
tion statistics of detectable objects (Watters et al. 2009). 

Also, the predicted pulse profile complexities for the two 
models differ. OG models are complex for a ~ £ ~ 90°; TPC 



models have complex profiles at intermediate a. In general, 
the radio-faint pulsars are double-peaked in both models, but 
the TPC picture has relatively more single- and zero-peaked 
solutions. 

One area of agreement between the TPC and OG models is 
the preponderance of large /n; this is a principal result of our 
study: inferred efficiencies should be increased by roughly 
an order of magnitude from those commonly assumed. In 
practice this means that some young pulsars must have high 
(> 10%) 7-ray efficiencies. This argues that a study of the 7- 
ray emission is an even better probe of magnetosphere struc- 
ture and energetics than previously thought. Our /□ correc- 
tions tighten the well-known trend for efficiency to scale as 
77 oc E^ 1 / 2 , although a few outliers remain. If, as more pul- 
sars are discovered by the LAT, we find that this trend re- 
mains strong, we may even be able to use observed fluxes, 
after /$i correction, for "standard candle" estimates of pulsar 
distances. Such arguments suggest relatively low distances 
for PSRs J2021+3651 and B1706-44. 

Importantly, for models with w(E) increasing with age, 
fn tapers off at small E. Thus the sky coverage decreases 
smoothly into a "death zone" where the gap shuts off. Near 
shut-off, the small solid angle on the sky ensures that old pul- 
sars appear relatively bright and are visible to large d. This 
will be an important signature in the population sums (Watters 
et al. 2009). Examination of the distribution of pulse proper- 
ties over the simulated pulsar sample is thus the key to testing 
the models and to producing the corrected luminosities that 
will make LAT data powerful in probing pulsar populations 
and spindown. 
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APPENDIX 



Here we show a selection of light curves for the two-pole caustic (TPC) and outer gap (OG) models. For both, we arrange the 
sample light curves on an a, (, grid so that the results may be compared with the "atlas" figures (Figures |2]|4|. Light curves are 
shown for four w values (0.01, 0.1, 0.2, 0.3). The number of peaks (total and broad only) and the maximum peak separation A 
(expressed here in %) as flagged by the automatic peak finder are given near the individual light curves. 

These light curves allow the reader to visualize the topological parameters summarized in the color "atlas" figures. For ex- 
ample, in the TPC plot one can see the broad, sinusoidal profiles that arise at low a and £ and thus earn the "zero major peak" 
classification. One also sees how the TPC model has complex pulses near a ~ 40°, £ ~ 60° and that these sometimes give small 
separations between the strongest peaks. For the OG plots, one sees that for a > 60° and small £ the normal first peak is often 
lost; this is a result of the rapid spreading of the high-altitude field lines contributing to this peak, which prevents formation of 
the sharp leading caustic. Also, one sees that the pulse shapes are often complex near £ = 90° and that the pulse shapes simplify 
toward basic double pulses as w becomes large. 
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FIG. 9. — Collection of sample light curves for the TPC model. Four select w (values in bottom right panel) are shown for each panel; the radio pole has closest 
approach at phase=0. The values for the number of all peaks, the number of broad peaks and the maximum peak separation (in %) are indicated by each curve. 
Intensities are normalized to pulse maximum. 
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FIG. 10. — OG model light curves, as for Figure[9] Four select w (values in bottom right panel) are shown for each panel; the radio pole has closest approach 
at phase=0. Intensities are normalized to pulse maximum. For small a or f , especially at large to, no pulses form. 



